Effect of intersubsystem coupling on the geometric phase in a bipartite system 
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The influence of intersubsystem coupling on the cyclic adiabatic geometric phase in bipartite 
systems is investigated. We examine the geometric phase effects for two uniaxially coupled spin— | 
particles, both driven by a slowly rotating magnetic field. It is demonstrated that the relation 
between the geometric phase and the solid angle enclosed by the magnetic field is broken by the 
spin-spin coupling, in particular leading to a quenching effect on the geometric phase in the strong 
coupling limit. 
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The geometric phase, originally conceived by Berry P| 
for cyclic adiabatic evolution of pure quantal states, has 
been extensively studied and generalized, for exam- 
ple to nonadiabatic evolution 0] , mixed states 0] , and 
open systems . The appearance of the geometric phase 
in composite entangled systems with no intersubsystem 
coupling has been analyzed for a pair of entangled spins 
in a time-independent uniform magnetic field '8| and for 
the case of entangled spin pairs in a rotating magnetic 
field i; it has also attracted interest in connection to 
polarization-entangled photon pairs ) the topology of 
the SO (3) rotation group Bell's theorem ^21 j ^ "^cU 
as in relation to the mixed state geometric phases of the 
subsystems [T^ . 

The importance of geometric phases to robust control 
of quantal systems, such as in fault tolerant quantum 
computation Q , has triggered extension of the geomet- 
ric phase for composite systems to interacting subsys- 
tems. This topic has been addressed for two isolated 
interacting spins (3, for systems entangled with a quan- 
tized driving field p^ . in relation to unitary representa- 
tions of quantum channels IT^ . and more recently for a 
pair of interacting spin— i particles with one of the spins 
driven by a slowly rotating magnetic field Concern 
about the effect of interaction on the geometric phase 
may also arise in the application of the geometric phase to 
systems with intra- variable couplings, such as, e.g., spin- 
orbit coupling in atomic systems, where the entanglement 
among a distinguished set of observables becomes an at- 
tractive issue in quantum information processing p^ . 

In this work, we develop the theory geometric phases 
of bipartite systems with intersubsystem coupling under- 
going adiabatic cyclic evolution. We examine the effect 
of intersubsystem coupling on the pure state geometric 
phase of the composite system, as well as on the geomet- 
ric phases associated with the reduced density operators 
of the subsystems. We calculate and analyze the geomet- 
ric phases in the case of two uniaxially interacting spin-^ 
(qubit) systems with the same magnetic dipole moment 
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and driven by a slowly precessing magnetic field; a case 
of relevance to, e.g., entanglement creation by adiabatic 
passage techniques [l9j as well as to NMR quantum com- 
putation [23, 1^ . Finally, we briefly discuss a possible 
extension of this analysis to systems with intra-variable 
coupling. 

Let a quantal system S be exposed to the Hamiltonian 
H(Q), Q being some external control parameters. Sup- 
pose that the Hilbert space TC of S is N dimensional and 
that Q varies around a closed path C : t € [0, T] — > Qt 
in parameter space, so that H{Qt) — H{Qo)- Expan- 
sion of the solution of the Schrodinger equation in the 
instantaneous eigenstates \n{Qt)) of H{Qt) yields 

N 
n=l 

If T is large enough, the adiabatic theorem |22| entails 
that transitions between the instantaneous energy eigen- 
states are negligible, making \cn\ approximately time- 
independent, so that the final state reads (h = 1 from 
now on) 

N J. 

\^{T)) = J2 Ue-Jo ^"(*)^*e^^"['^V(Qo)), (2) 

n=l 

where En (t) and 7„ [C] are the instantaneous nondegener- 
ate energy eigenvalue and the cyclic adiabatic geometric 
phase, respectively, associated with the nth energy eigen- 
state. 

Now, let us focus on the case where S has a natural 
bipartite decomposition in terms of subsystems Sa and 
Sb with corresponding Hilbert spaces Ha and Hb- For 
such S, the energy eigenvectors may be put on Schmidt 
form 



HQ,9)) - E '\Jpi"\Q,g)\ai"\Q,g)) <E> \bi"\Q,g)), 

k=l 

(3) 

where the Schmidt vectors ja^,"-* (Q, 17)) (Q, 5)) 
are characterized by {al^\Q , g)\a'f'\Q , g)) — Ski and 
{b'"^\Q , g)\b'l"'\Q , g)) — Ski, g is some set of fixed cou- 
pling parameters, and N = min {dimHa,dimnb). The 



2 



Schmidt decomposition is unique provided the nonvan- 
ishing coefRcients p)^^ are nondegenerate, i.e., p^"'' ^ 

(n) 

PI , Vfc, I. The geometric phase associated with the path 
C in parameter space may be written as 



7(^)[C;ff] = ln[C]=ij^dQ-{n{Q,g)\^Q\n{Q,g)) 

- E(rS[C;.9] + fiS[C;5]), (4) 



(n) 

where we have used that the p), 's sum up to unity. Here 



«;fc 



[C;g]=z(^dg.(Ci"^(Q,g)|VQCi"^(Q,.g)) 



(5) 



are geometric phases of the weighted Schmidt vectors 

\^u\Q^9)) = [p'':\Q,9)f\t\Q^9)) pertaining to 
subsystem S^=a,b- 

Next, let us introduce the concept of nontransition 
eigenstates. These are defined as energy eigenstates 
where the p^'s are time-independent. For vanishing in- 
tersubsystem couphng, i.e., when g — 0, only such states 
occur since the time evolution operator then takes the 
bi- local form Uab = Ua'S'Ub, which exactly preserves the 
Schmidt coefRcients p^"^ On the other hand, transitions 
usually occur in the presence of coupling and the non- 
transition condition is in general only valid for specific 
paths. Closed paths of this kind are rare but have been 
found and studied for a spin-spin interaction model in 
Ref. [13. 

Suppose that there exists a nontransition state tracing 
out a closed path V in parameter space. For such a path, 
we may compute the cyclic geometric phase as 

ji-^V-g] ^ ^Y.Pk\Qo.9)(T^:l['D;g]+Ti-hV;g] 



(6) 



where 



rg[I?;g] = ^j^dQ•(ei"^Q,5)|VQCf^(Q,5)) 

= [pi"^(Qo,5)]"'rg[P;5] (7) 

with ^ — a,b, constitute the one-particle geometric 
phases of the Schmidt vectors. The nontransition prop- 
erty makes it natural to extend Ref. @ and define mixed 
state geometric phases for the two interacting subsystems 



7^ ' [V; g] = arg ^ pf' (Qo, s) exp [V; g]) . (8) 

k 

Thus, the geometric phases of the subsystems are taken 
as the average of phase factors pertaining to the non- 
transition eigenstates of the reduced density operators, 
weighted by the corresponding eigenvalues. 

Now, consider two qubits a and b as represented by a 
pair of spin-i particles coupled by a uniaxial exchange 



interaction in the z direction. In the presence of a 
time-dependent external magnetic field B(i) — Boh{t) 
with the unit vector n = (sin 9 cos 6, s in 9 sin (j>, cos 6), the 
Hamiltonian of this system reads |23| 



Hit) - AJS'^S^ + pBit) ■ [Sa + Sfc 



(9) 



with fj, the magnetic dipole moment assumed to be equal 
for the two spins. The first part of the Hamiltonian de- 
scribes the exchange interaction (spin-spin coupling) with 
coupling constant J, fj, is the gyromagnetic ratio, and 
= {S^,S^, S^) is the $th spin operator (^ — a, b). An 
explicit physical scenario for this model could be NMR 
experiments on -'^^C-labeled trichloroethylene in which 
the nuclear spins of the two ^^C nuclei could act as two 
spin— i systems with nearly the same magnetic dipole 
moment (the chemical shift of the two nuclei for this sys- 
tem is typically a fraction 10~^ of their precession fre- 
quency j2l|). In terms of the orthonormalized total spin 
eigenstates \S; M), 5* = 0, 1 and M = -S, . . . , 5, in the 
z direction and in units of p,Bo , the Hamiltonian Eq. 
can be expressed in the block-matrix form 



Hit) 



Hcit) 








(10) 



where 




^ sin 6*6"^ 







-9 TTs sm 



^sin6le"*'^ 

v2 



V2- 

g + cos ( 



pi0 



(11) 

with rescaled coupling constant g — -jf^ that may be 
controlled by changing the magnitude of the external 
magnetic field. Thus, the spin singlet |0;0) is decoupled 
from the triplet states |1;M), M = -1,0, 1. Therefore, 
in contrast to the case where only one of the subsystems 
interacts with B(i) [l^, the adiabatic geometric phase 
acquired by the singlet state vanishes here. 

First, let us consider the case of vanishing spin-spin 
coupling characterized hy g = Q. In this case, the total 
spin in the n direction commutes with Hit) so that the 
adiabatic geometric phases may be expressed in terms of 
the total spin projection quantum number M along the 
magnetic field and the solid angle Vl enclosed by the path 
T) in parameter space. For M = ±1 we obtain 



lab 



(12) 



where the former follows from the standard Berry for- 
mula —MVt 1]. The ill = eigenstates are two-fold 
degenerate and = ajl; 0) -I- /3|0; 0) for any complex 

numbers a and (3 is an energy eigenstate. While 7^^-* [D; 0] 
vanishes when taking |\E'*^°)) around V, the correspond- 
ing mixed state geometric phase for the two subsystems 
become 



3 



= — arctan 



2Re(a*/3) tan 



n 



(13) 



provided 2Re(a*/3) ^ 0. On the other hand, when 
2Re(Q;*/3) = 0, which for example occurs if the two spins 
are associated with indistinguishable entities, in case of 
which the singlet and triplet states cannot mix, the re- 
duced density operators of the subsystems are degener- 
ate, and the corresponding geometric phases become un- 
defined since no direction in space is singled out by the 
corresponding Bloch vectors. In all other cases, we have 
7(f) [D; 0] = 7^*^^ [D; 0] + 7^^^^ [V; 0], which is due to the 
spherical symmetry of the model in the 5 = case. 

For (7^0, the spherical symmetry is broken, and 
there is no simple relation neither between the geomet- 
ric phases and the solid angle nor between the geometric 
phase of the composite system and those of the subsys- 
tems. Thus, to proceed we need to diagonahze Hc{t) to 
obtain its eigenstates as (n = — , 0, 



e^^yl(")(0,g)|l;-l>+i3(")(e,.g)|l;O) 
+e-^^C(")(0,5)|l;l) (14) 



with 

_B(«) 



1 



1 



- cos 61] sine, 



1 



[(X(«))2 

-f cos 61] sin 6*, 



V2Af(") 

(X(«))4 + (1 - 3cos2 6l)(X("))2 + cos^ 0.(15) 



The shifted instantaneous energy eigenvalues X^") = 
X(")(0,5) = En - g oi the Hamiltonian Eq. are 
solutions of 



2gX^ - X-2g cos^ 6* = 0, 



(16) 



which yields X^^\e,g) ^ ±cos9 and X^^\9,g) = -2g 
in the limit of g — > 00. The Schmidt coefficients read 



(") 
Pi 



1 (") ^ I 1 

= 1 - = 2 ( 1 



^2(B("))2 + (c(«) - 



(17) 

Thus, p^"' and determined by the (f) indepen- 

dent effective Bloch vector r^") and it follows that the 
nontransition paths are those where 6 is constant. For 
closed paths V : t e [0,T] {(t)t,Ot) = {2Trt/T,9) with 
j-in) Q J" large, we uniquely obtain the adiabatic 
geometric phases for the corresponding Schmidt vectors 
pertaining to subsystem 6*5=0, & as 



ri':l[D,g] 



= -rg [P, g] = -TT (1 - F(") cos 9) , (18) 




0.5 1 1.5 2 2.5 3 1 2 3 4 5 6 

Coupling constant g 

FIG. 1: An illustration of the geometric phases of the com- 
posite system (nc), the sum of the corresponding mixed state 
geometric phases for the two qubits (nm), and the geomet- 
ric phase pertaining to one of the Schmidt vectors (ns), cor- 
responding to the three instantaneous eigenvectors l^p'"^), 
where n — — ,0, + and |\I>'"') become the spin triplet states 
|1; —1), 11; 0), |1; 1) in the strong spin-spin coupling limit. The 
horizontal axes show the rescaled dimensionless coupling con- 
stant g. The composite system is assumed to undergo an 
adiabatic evolution in one of the instantaneous eigenstates, 
such that the corresponding Schmidt coefficients are time- 
independent. These states are characterized by constant po- 
lar angles (a) 9 = n/6, (b) 9 = it/3, (c) 9 — 7r/30 (close to 
zero) and (d) 9 = 97r/20 (close to 7r/2). 



where the scale factor 

pM = F^"'>{9,g) 



sm9 



V(X("))4 - 2(X("))2 cos2 9 + cos2 9 



(19) 



comprises the effect of intersubsystem coupling on the 
geometric phase of the Schmidt vectors. The geometric 
phases of the composite system and those of the subsys- 
tems read 



27rr(") f 1 - cos 6'), 



7^"^[X>,g] = -arctan r("Han (7r[l - i^^") cos6']) 

(20) 

with ^ = a,b. 

The dependence of the geometric phase upon the cou- 
pling constant g is illustrated in Fig. An interesting 
feature of Fig. ^ is that all geometric phases for the 
composite system and its subsystems tend to an integer 
multiple of 27r when 5 — > oo. This limit corresponds 
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FIG. 2: The geometric phase in units of tt for the eigenstates 
1^^*') and l^^"') vs. the rescaled dimensionless coupling con- 
stant g with 9 = TT sin(7rt/T), T ^ 1 being the period of 
precession chosen in accordance with the requirement of adi- 
abaticity. 

to the case when the second term in the Hamiltonian 
in Eq. ^ can be ignored. In other words, the triplet 
states |1; M) would become the instantaneous eigenstates 
of the system with g oo, thus making the enclosed 
area in state space to vanish and thereby the geomet- 
ric phase factors become trivial. It is worth stressing 
that the Hamiltonian in Eq. ^ has permutation sym- 
metry, resulting in 7^"^ = 7^"'[I?;(7]. In fact, the 
mixed state geometric phases presented in Fig. ^are for 
7^"^ [V; g]+jjr^ [V- g] = 2jt^ [V- g] = 27^"^ [V- g] ; the sum 
of mixed state geometric phases of the two qubits being 
equal to twice the geometric phase for each subsystem. 
Note in particular that 7^^^^ [V; g] ^ 7!"^ [V; g] -^7^"^ [V; g] 
in general. The state of the composite system are en- 
tangled at most time when g ^ 0, this indicates that 
there are at least two nonvanishing Schmidt coefficient 
Pk and thus most the geometric phases for the Schmidt 
vector in Eq. jT)) would have similar properties in the 
limit g 00 except that pertaining to j^''"'), where 
X(")(6',5 00) —g, and then F^"^ tends to zero, 

consequently, T^^^i'^Td ^ 00] ^ tt as shown by Os lines 
in Fig. H] 

For transition paths C along which the polar angle 9 
and thereby the Schmidt coefficients vary, the geometric 
phase of the composite system reads 

^ab[C,g] = -j^r(")(l-F(")cos^^)d0. (21) 



In Fig. 2, 7^" [C, g] is shown as a function of the rescaled 
coupling constant g for the closed path C : t G [0,T] ^ 
{(j)t,et) = (7rt/T,7rsin(7ri/r)). With g ^ 00, the phases 
tend to zero, in analogy with the nontransition case dis- 
cussed above. 

An interesting extension of above analysis is to the case 
of systems with intra-variable coupling. For example, 
one may consider an atom with electronic orbital and 
spin angular momentum L and S, respectively, precessing 
in a time-dependent magnetic field B, the Hamiltonian 
describing such a system reads, H = iin{t) ■ (L + 2S) + 
(/L-S, the last term describes the spin orbit coupling. An 
analysis is expected to show that the spin-orbit coupling 
would affect the geometric phase of the atom in a similar 
way as in the intersubsystem coupling case. 

In conclusion, we have analyzed the cyclic adiabatic 
geometric phase of bipartite systems, focusing on the ef- 
fect of intersubsystem coupling. We have distinguished 
two different kind of evolution in regard to whether or 
not the Schmidt coefficients are time-dependent. The 
geometric phases of the subsystems naturally extends in 
terms of the standard mixed state geometric phase in 
the nontransition case, i.e., when the Schmidt coefficients 
and thereby the eigenvalues of the corresponding reduced 
density operators are fixed. We have found a striking ev- 
idence for a strong influence of the intersubsystem cou- 
pling on the geometric phases for two uniaxially coupled 
spin— i systems, such as a divergence from the standard 
relation to the solid angle enclosed by the driving mag- 
netic field, leading to a quenching effect on the geometric 
phases in the strong coupling limit. This latter result has 
also been demonstrated in the transition case, where the 
Schmidt coefficients are changing around the curve in pa- 
rameter space. Physically, the quenching effect may be 
viewed as a consequence of the broken spherical symme- 
try as expressed by the existence of a preferred direction 
in space singled out by the uniaxial exchange direction 
of the spins, making the energy eigenstates independent 
of the slowly rotating magnetic field. 
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